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1 .  INTRODUCTION 


In  a  series  of  papers  [1]-[6],  researchers  at  the  NOAA  Wave  Propagation 
Laboratory  have  described  the  operation  of  a  dual  frequency  (20.6  and  31.6 
GHz)  radiometer  and  the  procedure  for  retrieving  the  precipitable  water  vapor 
(V)  and  liquid  (L)  in  a  non-precipitating  atmosphere  using  this  radiometer. 
Statistical  retrieval  algorithms  are  used  for  the  determination  of  V  and  L. 

In  ♦'his  paper  we  describe  a  method  which  can  give  estimates  of  V  and  L  in 
a  non-precipitating  atmosphere,  and  for  V  in  light  rain,  from  emission  meas¬ 
urements  using  a  three-channel  radiometer  operating  near  the  22.235  GHz  water 
vapor  line. 


E 


2.  PRINCIPLES  OF  OPERATION 


The  attenuation  in  a  non-precipitating  atmosphere  at  the  radio  wave¬ 
lengths  of  interest  in  this  study  is  attributable  to  absorption  by  water 
vapor,  oxygen,  and  the  water  droplets  in.  clouds  and  fog.  The  absorption 
coefficients  for  water  vapor  and  oxygen  are  found  in  the  literature,  e.g., 
Waters  [7]  and  Liebe  [8].  For  water  vapor, 

-644/T 

C  0 

“w  “  C1p  Avw  ^  “ — j -  f(v,Avw)  +  c_]  km  1  (1) 

where  p  is  the  water  vapor  density  in  gm~3,  v  is  the  frequency  in  GHz,  Av 

w 

is  the  linewidth  in  GHz,  v  =  22.235,  T  is  the  temperature  in  K,  and  f(v,Av  ) 

o  w 

is  the  Van  Vleck-Weisskopf  line  shape.  For  the  oxygen  molecule, 

“o*  -  i.n  p  r3  v  I  [|UNJ2  v„,  r  <»,v„,>  ♦ 
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F(  v )  ] 


"EN/kT 


km 
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where  P  is  the  total  pressure  in  mb,  |uN+|2  are  the  matrix  elements  for  the 
AJ  =  ±  1  transitions,  |uNo|2  are  the  matrix  elements  for  the  non-resonant 
absorption,  vN+  are  the  resonant  frequencies  for  these  transitions,  f(v,vN+) 
are  the  corresponding  kinetic  line  shapes,  F(v)  is  the  non-resonant  line 
shape,  and  is  the  energy  for  the  rotational  state.  The  absorption 
coefficient  for  oxygen  is  obtained  by  summing  over  all  transitions  which  have 
contributions  at  the  frequencies  of  interest.  The  cloud  or  fog  absorption  [9] 
is  given  by 

a  =  1.885  y  Im  (-  ^ 0  km_1  (3) 

c  A  m2  +2 

where  M  is  the  liquid  water  content  in  gm~3,  \  is  the  wavelength  in  cm,  m  is 
the  complex  index  of  refraction  of  liquid  water,  and  Im  means  the  imaginary 
part.  The  index  of  refraction,  m,  depends  on  the  temperature  and  wavelength 
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such  that  the  cloud  absorption  coefficient  is  very  nearly  proportional  to 

[10]. 


The  rain  attenuation  coefficient  is  given  by 


/  Q  (D)  N  (D)  dD 
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where  D  is  the  drop  diameter,  Qt(D)  is  the  extinction  expressed  in  terms  of 

the  Mie  coefficients,  and  N(D)  is  the  drop  size  distribution.  For  a  given 

x 

temperature  and  rain  rate,  o^,  -  v  .  Near  the  water  vapor  line,  x  varies  in 
the  range  1.9  -  2.6  for  temperatures  between  0  and  10  degrees  Centigrade  and 
rain  rates  between  0.25  and  10  mm/hr  for  the  Laws  and  Parsons  [11],  and  the 
Marshall  and  Palmer  [12]  drop  size  distributions. 

Estimates  of  the  precipitable  water  vapor  and  liquid  can  be  obtained  by 
the  following  procedure.  First  we  assume  that  at  a  given  elevation  angle  the 
total  opacity,  t,  can  be  determined  by  emission  measurements.  In  this  study 
the  attenuation  direction  is  assumed  to  be  the  zenith.  Opacities, 
t i  ,  t2,  are  determine';1  at  three  frequencies,  Vl  ,  v2,  and  v  where 

v2  is  the  line  center  of  the  water  vapor  line  at  22.235  GHz.  The 
frequencies  v1  and  are  in  the  vicinity  of  and  are  chosen  such  that 


21,2  2, 

=  -  (v1  +  ) 


(5) 


We  show  that,  for  practical  purposes,  the  derived  quantity, 


Td  =  t2  "  2  (T1  +  t3} 


(6) 


is  dependent  only  on  some  integrated  function  of  the  water  vapor  resonant 
absorption  at  22.235  GHz.  The  utility  of  the  method  depends  on  whether  is 
a  good  measure  of  V. 

The  opacities  t2>  and  are  estimated  from  the  radiometric  formula, 
given  by 
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T 

t  =  In  — —  (7) 

m  e 

where  Tm  is  the  effective  medium  temperature,  and  Te  is  the  emission  tempera¬ 
ture  of  the  atmosphere.  In  practice  Te  would  be  the  measured  emission  temper¬ 
ature,  and  Tm  would  be  an  empirically  derived  constant,  dependent  primarily  on 
local  meteorological  conditions  and  the  characteristics  of  the  antenna  sys¬ 
tem.  The  quantity  Tm  would  be  statistically  determined  from  emission  and 
attenuation  measurement  on  cloudless  days  using  an  extra-terrestrial  source  or 
by  comparing  with  a  previously  calibrated  instrument.  At  the  same  time, 
would  be  derived,  and  correlated  with  the  precipi table  water  vapor,  V, 
which  would  be  determined  from  radiosondes  or  other  calibrated  instruments. 
Also,  V  would  be  correlated  with  the  opacities  ,  x^,  and  x^*  From  these 
measurements  on  cloudless  days,  statistical  relations  between  x^  and  V,  and  V 
and  xQ  (clear  atmosphere  opacity)  are  obtained.  In  operation,  the  cloud 
opacity  is  obtained  from 


and  the  precipitable  liquid  is  obtained  from  tc  using  estimates  of  the  cloud 
temperature . 

In  the  operation  given  by  (6)  any  opacity  with  a  monoto.nic  behavior  in 
the  frequency  range  v1  to  would  tend  to  cancel  (exactly  for  a  v  or  v  de¬ 

pendency).  In  fact,  the  only  kind  of  function  which  would  not  cancel  is  one 
with  a  turning  point,  such  as  a  resonant  absorption;  also,  because  (6)  is  a 
differential  operation,  the  error  in  x^  is  relatively  insensitive  to  errors  in 
Tm,  and  inferentially ,  relatively  insensitive  to  the  temperature  distribu¬ 
tion.  Near  the  water  vapor  line  the  oxygen  absorption  coefficient  is  approx¬ 
imately  a  linear  function  of  frequency  and  is  sufficiently  small,  that  in  (6) 
for  the  five  model  atmospheres  to  be  used  as  examples,  the  oxygen  opacities 
are  cancelled  to  a  value  less  than  10-^.  Near  the  water  vapor  line,  for 
temperatures  -10  to  +10  degrees  Centigrade,  the  cloud  absorption  coefficient 
varies  as  vX  where  x  ranges  from  1.88  to  1.95.  For  a  1  km  thick  cloud  with  M 
=  1  gm-3  in  the  above  temperature  range,  the  cloud  opacities  in  (6)  cancel  to 
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a  value  less  than  10-^.  For  a  3  km  thick  rain  layer  with  temperatures  in  the 
range  0  to  10°C,  and  R  in  the  range  0.25  to  10  mm/hr,  the  rain  opacities 
cancel  to  a  value  less  than  5  x  10”^  for  the  two  drop  size  distributions 
mentioned  previously. 


3.  SAMPLE  APPLICATION 


The  procedure  outlined  in  the  previous  section  will  be  used  to  retrieve  V 
and  L  in  five  supplemental  atmospheres  given  in  the  Handbook  of  Geophysics  and 
Space  Environments  [13]  in  a  simulated  operation.  These  atmospheric  models 
are  the  (1)  tropical;  (2)  subtropical,  Jul .  ;  (3)  subtropical,  Jan.;  ( -4 ) 
midlatitude,  Jul.;  and  (5)  midlatitude,  Jan.  The  top  of  these  atmospheres  was 
assumed  to  be  21  km,  and  the  relative  humidity  above  10  km  was  taken  to  be  the 
same  as  that  at  1 0  km  given  in  the  Handbook.  A  1  km  thick  cloud  layer,  with  M 
=  0.25,  0.50,  and  1.0  gm-^  with  the  cloud  bottom  at  2  or  3  km  was 
superimposed  on  these  models. 

Emission  temperatures  and  opacities  were  calculated  for  the  five  atmos¬ 
pheric  models  at  three  frequencies,  v1  =  20  GHz,  =  22.235  GHz,  and  v,  = 
24.265  GHz.  The  choices  of  v1  and  are  not  critical,  as  long  as  (5)  is 
satisfied  and  the  frequencies  are  spaced  at  a  reasonable  fraction  of  a  line- 
width.  These  calculations  were  made  with  and  without  cloud  cover.  The  total 
opacities  t^,  and  were  derived  from  (7)  assuming  that  the  calculated 
emission  temperature  Te  to  be  the  measured  one.  In  this  study  we  assume  the 
effective  medium  temperature  Tm  to  be 


T  -  T  -  13  (9) 

m  o 

where  T0  is  the  surface  temperature.  In  effect  this  assumes  an  infinitely 
narrow  antenna  beamwidth  with  100?  antenna  beam  efficiency.  However,  these 
assumptions  do  not  affect  the  principles  of  operation  and  the  conclusions 
reached  about  the  retrieval  method.  Using  the  value  of  Tm  in  (9),  the  differ¬ 
ence  between  the  derived  opacity  using  (7)  and  the  calculated  opacity  was  less 
than  0.005  (0.02  dB)  for  all  models,  with  and  without  cloud  cover. 

In  retrieving  V  and  the  clear  atmosphere  opacity,  linear  relations  of  the 

form 


V 


Vd 


(10) 


9 


and 


Ks 


L01 


b0  +  b!  v 


(ID 


are  assumed.  The  coefficients  a-p  bg,  and  b^  were  determined  from  the  best 


fit  linear  regression  analysis  using  the  calculated  values  of  t 1 ,  i^t  and  V 


in  the  five  cited  model  atmospheres  without  clouds.  In  practice,  Tm,  a1 ,  b0, 
and  b-|  would  be  statistically  determined  by  correlating  emission  measurements 
at  a  single  site  on  cloudless  days  with  independently  measured  values  of  V  and 
Tq.  In  (11)  Tqi  is  the  estimated  clear  atmosphere  opacity  at  ^  .  Any  of  the 
other  frequencies  could  have  been  utilized. 

~3 


Results  are  shown  in  Figs.  1  and  2  for  a  1  km  thick  cloud  with  M  =  1  gm 
superimposed  on  the  five  atmospheric  models.  The  cloud  bottom  was  at  3  km. 
Similar  results  were  obtained  with  the  other  cloud  conditions.  The  triangles 


are  the  calculated  values  of  V  and  (Fig.  1),  and  tq1  and  V  (Fig.  2)  for  the 


five  numbered  cloud-free  model  atmospheres  listed  in  the  beginning  of  the  sec¬ 
tion.  In  Fig.  1  the  circles  give  the  retrieved  values  of  V  using  the  proce¬ 
dure  described  earlier.  The  y-difference  between  the  triangles  and  circles  is 
the  error  in  the  retrieved  values  of  V.  The  largest  error  occurs  for  model  1 , 
which  had  a  water  vapor  distribution  shape  that  was  significantly  different 
from  the  others.  In  this  model,  there  was  relatively  more  water  vapor  lower 
in  the  atmosphere.  This  can  be  seen  in  Fig.  3.  Because  the  water  vapor  line 
is  pressure  broadened,  is  sensitive  to  the  distribution  of  the  water  vapor 
density.  For  example,  it  is  seen  that  if  the  water  vapor  is  distributed  at  a 


higher  altitude,  less  V  is  required  for  a  given  than  for  the  case  where  the 


water  vapor  is  distributed  lower  in  the  atmosphere.  In  Fig.  2  the  circles 
indicate  the  retrieved  value  of  the  cloud-free  atmosphere  opacity,  using  the 
retrieved  value  of  V  in  (11).  The  y-difference  between  the  triangles  and 
circles  is  the  error  in  the  estimated  cloud  opacity  for  that  model. 


At  frequencies  near  the  water  vapor  line  the  attenuation  through  light 
rain  is  primarily  absorptive  and  the  effective  medium  temperature  used  in  the 
radiometric  formula  for  non-raining  conditions  can  be  used  to  determine  the 


total  atmospheric  opacity  in  light  rain.  Also,  as  mentioned  before,  is  a 


differential  quantity  and  errors  in  are  relatively  insensitive  to  small 


errors  in  T_.  Using  this  procedure  in  light  rain  will  give  estimates  of  V. 
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Relation  Between  Preci pi  table  Water  Vapor  and  the  Derived 
Differential  Opacity  where  the  Numbers  Indicate  the  Atmospheric 
Models.  The  triangles  are  calculated  and  the  circles  are  the 
retrieved  values 


Relation  Between  the  Clear  Atmosphere  Opacity  and  Precipitable 
Water  Vapor  for  the  Indicated  Atmospheric  Models.  The 
triangles  are  calculated  and  the  circles  are  the  retrieved 
values . 


H.  RESULTS  AND  CONCLUSIONS 


The  accuracy  of  this  method  for  retrieving  V  and  L  depends  on  the  follow¬ 
ing  assumptions:  (1)  the  total  opacity  can  be  determined  by  emission  measure¬ 
ments;  (2)  the  differential  opacity,  gives  a  measure  of  the  precipitable 
water  vapor;  and  (3)  the  precipitable  water  vapor  is  a  measure  of  the  clear 
atmosphere  opacity  in  non-raining  conditions.  Assumptions  (1)  and  (3)  are 
essentially  utilized  in  the  dual  radiometer  system  operating  at  20.6  and 
31.6  GHz.  It  was  shown  that  is  a  function  only  of  the  water  vapor 
absorption  and  is  sensitive  to  the  height  distribution  of  the  water  vapor. 

If,  at  a  given  site,  the  form  of  the  water  vapor  distribution  is  highly 
variable,  the  rms  errors  in  the  retrieved  V's  will  be  correspondingly 
increased.  In  this  case,  instead  of  a  simple  t^-V  relation,  more  channels  can 
be  utilized  in  the  v1  -  interval  and  a  retrieval  algorithm  based  on 

a  t  .  curve  can  be  imolemented.  The  Tj  curve  also  gives  some  information  about 
a  d 

the  water  vapor  distribution. 

It  is  not  possible  to  give  a  quantitative  assessment  of  the  accuracy  of 
this  method  for  retrieving  V  and  L,  in  the  limited  number  of  illustrated 
examples.  However,  the  method  does  have  several  features  which  make  it 
attractive  for  field  use.  It  utilizes  a  single  radiometer  and  its  output  is  a 
measured  quantity  which  is  a  function  only  of  the  water  vapor  density.  The 
system  can  retrieve  V  in  light  rain,  and  with  additional  channels,  can  provide 
information  about  the  wate^  vapor  distribution. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff’s  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  therroa1  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exciraer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of -f ield-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards ,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  faul t -tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communlcat ion  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  microwave/millimeter  wave  measurements, 
diagnostics  and  radiometry,  microwave/millimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnet ic 
propagation  phenomena,  space  commun icat ion  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  an<T  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy -induced 
env i ronments . 

Space  Sciences  Laboratory:  Magnetospher ic,  auroral  and  cosmic  ray 
physics,  wave -part ic le  interactions,  magnetosphe r ic  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomv, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumental  ion. 
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